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, Abstract 

In this paper, we study the existence of the density associated to the exponential 
functional of the Levy process £, 



•— » I eq := / e«' da, 

where e q is an independent exponential r.v. with parameter q > 0. In the case when 
£ is the negative of a subordinator, we prove that the density of I eq , here denoted 
by k, satisfies an integral equation that generalizes the one found by Carmona et al. 
,-h [7]. Finally when q = 0, we describe explicitly the asymptotic behaviour at of the 

density k when £ is the negative of a subordinator and at oo when £ is a spectrally 
positive Levy process that drifts to +oo. 
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O 1 Introduction 

A real-valued Levy process is a stochastic process issued from the origin with stationary 
and independent increments and almost sure right continuous paths with left-limits. We 
write £ = (£t,t > 0) for its trajectory and P for its law. The law P of a Levy processes 
is characterized by its one-time transition probabilities. In particular there always exists 
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a triple (a, a 2 , II) where a G K, a 2 > and II is a measure on M\{0} satisfying the 
integrability condition L(l A x 2 )II(dx) < oo, such that, for £ > and zel 

E[e iz ^}=exp{-^(z)t}, (1.1) 

where 



Sl/fz) = iaz 



+ -<rV + / (l-e^ + ^xl { | x | <1} )n(dx). 



In the case when £ is a subordinator, the Levy measure II has support on [0, <x>) and 
fulfils the extra condition J, Q ,(1 A a;)II(dx) < oo. Hence, the characteristic exponent ^ 
can be expressed as 

#(*) = -«C2 + / (l - e^)n(dx), 

</(0,oo) ^ ' 

where c > and is known as the drift coefficient. It is well-known that the function \I/ 
can be extended analytically on the complex upper half-plane, so the Laplace exponent 
of £ is given by 

0(A) := -logEfe-^ 1 ] = *(«A) = c\+ f (l - e- Ax )n(dx). 

J(0,oo) ^ ' 

Similarly, in the case when £ is a spectrally negative Levy process (i.e. has no positive 
jumps), the Levy measure II has support on (— oo,0) and the characteristic exponent \1/ 
can be written as 

W(z) = iaz + -a 2 z 2 + (l - e izx + izxl {x> _ 1} )u(dx). 

It is also well-known that the function \1/ can be extended analytically on the complex 
lower half-plane, so its Laplace exponent satisfies 

i){\) : = logE[e A?1 ] = -*(-iA) = a\+ \a 2 \ 2 + I (e Xx - 1 + -Xxl {x> . 1} )u(dx). 

2 J (-oo,0) v ' ' 

In this article, one of our aims is to study the existence of the density associated to the 
exponential functional 

I eq := / e ?s ds, 
Jo 

where e q is an exponential r.v. independent of the Levy process £ with parameter q > 

If q — 0, then e q is understood as oo. In this case, we assume that the process £ drifts 

towards — oo since it is a necessary and sufficient condition for the almost sure finiteness 

of I := /oo, see for instance Theorem 1 in Bertoin and Yor [3]. 

Up to our knowledge nothing is known about the existence of the density of I e when 

q > 0. In the case when q = 0, the existence of the density of I has been proved by 

Carmona et al. [7j for Levy processes whose jump structure is of finite variation and 

recently by Bertoin et al. [2] (see Theorem 3.9) for any real-valued Levy process. In 



particular when £ is the negative of a subordinator such that E[|£i|] < oo, Carmona et al. 
[7] (see Proposition 2.1) proved that the r.v. / has a density, k, that is the unique (up to 
a multiplicative constant) L 1 positive solution to the equation 

POO 

(l-cx)k(x) = Tl(\og(y/x))k(y)dy, xG(0,l/c), (1.2) 

J X 

where II (x) := II(a;, oo). 

Here, we generalize the above equation. Indeed, we establish an integral equation for 
the density of I e , q > 0, when £ is the negative of a subordinator and we note that when 
q — 0, the condition E[|£i|] < oo is not essential for the existence of its density and the 



validity of (1.2). 



Another interesting problem is determining the behaviour of the density of the expo- 
nential functional J at and at oo. This problem has been recently studied by Kuznetzov 
[13] for Levy processes with rational Laplace exponent (at and at oo), by Kuznetsov 
and Pardo [IE] for hypergeometric Levy processes (at and at oo) and by Patie [19] for 
spectrally negative Levy processes (at oo). In most of the applications, it is enough to have 
estimates of the tail behaviour P(J < t) when t goes to and/or P(J > t) when t goes to 
oo. The tail behaviour P(I < t) was studied by Pardo [20] in the case where the underlying 
Levy process is spectrally positive and its Laplace exponent is regularly varying at infinity 
with index 7 e (1, 2), and by Caballero and Rivero [6] in the case when £ is the negative 
of a subordinator whose Laplace exponent is regularly varying at 0. Furthermore, the tail 
behaviour P(J > t) has been studied in a general setting, see [8j [18j |21j [22] . The second 
main result of this paper is related to this problem. Namely, we describe explicitly the 
asymptotic behaviour at of the density of I when £ is a subordinator which in particular 
implies the behaviour of P(J < t) near 0. 

The paper is organized as follows: in Section 2 we state our main results, in particular 
we study the density of I e and the asymptotic behaviour at of the density of the 
exponential functional associated to the negative of a subordinator. Section 3 is devoted 
to the proof of the main results and in Section 4, we give some examples and some 
numerical results for the density of I eq when the driving process is the negative of a 
subordinator. 

2 Main results 

Our first main result states that I e has a density, for q > 0. Before we establish our 
first Theorem, we need to introduce some notation and recall some facts about positive 
self-similar Markov processes (pssMp) which will be our main tool in this first part. 

Let (£j, t > 0) be the process obtained by killing £ at an independent exponential time 
of parameter q > 0, here denoted by e q . The law and the lifetime of ^ are denoted by P^ 
and (3, respectively. 



We first note that 

-0 



(l,FA = ( I exp{& t }dt,P t J = ( / V'dt.PJ . 

For x > let Q x be the law of X^ x \ the positive self-similar Markov process with self- 
similarity index 1 issued from x, associated to £t via its Lamperti's representation (see 
[T7] for more details on this representation), that is for x > 

. ,. , , j x exp [C( t /x) } > if T ( t / X ) < °° 



where 



if r(t/x) = oo, 



r(s) = inf < r > : / e 5t dt > s L inf {0} = oo, 



and is understood as a cemetery state. The process X^ is a strong Markov process and 
it fulfills the scaling property, i.e. for k > 0, 

(kx<jlt>o)±(x?*\t>o). 

We denote by T (x) := inf{£ > : x\ x) = 0}, the first hitting time of X& at 0. Observe 
that for s > 0, we have the following equivalences, 

r 13 t 

r(s)<oo iff t(s)<(3 iff s< e^dt. 

Jo 

Hence, it follows from the construction of X that the following equality in law follows 

(t ,qi) = n\^dt,A. 

Now, we have all the elements to establish our first main result. It concerns the exis- 
tence of the density of I e . 

Theorem 2.1. Let q > 0, then the function 

1 



hit) 



X-Hteib} 



t > 0, 



is a density for the law of I e 



Corollary 2.2. Assume q > and that £ is a subordinator. Then the law of the r.v. I e 
is a mixture of exponential, that is its law has a density h on (0, oo) which is completely 
monotone. Furthermore, lirn^o h(t) = q. 



In the sequel, we assume that £ = —( where ( is a subordinator and we denote by 
U q (dx) the renewal measure of the killed subordinator (Q,t < e q ), i.e. 



E 



f(Ct)dt 



[0,00) 



f(x)U q (dx) 



(2.3) 



where / is a positive measurable function. If the renewal measure is absolutely continuous 
with respect to the Lebesgue measure, the function u q (x) = U q (dx)/dx, is usually called 
the renewal density. If q = 0, we denote Uq and w by U and u. 

Before stating our first main result, which is a generalization of the integral equation 



(1.2) of Carmona et al. for subordinators, in the next proposition we establish that the 
density of I e solves an integral equation in terms of the renewal measure U q . 

Proposition 2.3. Let q > 0. The random variable I e has a density that we denote by k, 
and it solves the equation 



k(x)dx 



k(ye x )U q (dx) 



almost everywhere. 



(2.4) 



The next result generalizes (1.2) 



Theorem 2.4. Let q > 0. The random variable I eq has a density that we denote by k, 
and it solves 



(l-cx)k(x)= Ti(\og(y/x))k(y)dy + q k(y)dy «6(0,l/c), (2.5) 

J X J X 

Conversely, if a density on (0, 1/c) satisfies this equation then it is the density of I e . 

The importance of the above result will be illustrated in Theorem 2.5 where we study 
the asymptotic behaviour at of the density k, and in Section 4 where we provide some 
examples where k can be computed explicitly. Further applications have been provided 
in Haas [H] and Haas and Rivero [12] where this equation has been used to estimate the 
right tail behavior of the law of I and to study the maximum domain of attraction of /. 



The following corollary is another important application of equation (2.5). In parti 



cular, it says that if we know the density of the exponential functional of the negative 
of a subordinator, say k, then for p > 0, x p k(x) adequately normalized is the density of 
the exponential functional associated to the negative of a new subordinator. The proof of 
this fact follows easily by multiplying in both sides of equation (2.5) by x p . Such result 



also appears in Chazal et al. [9] but in terms of the distribution of I e not in terms of its 
density. 

Corollary 2.5. Let q > 0, p > 0, c p a positive constant satisfying 

f 

i(0,oo) 



and suppose that when q > the renewal measure U q has a density. Then the function 
h(x) := c~ 1 x p k(x) is the density of the exponential functional of the negative of a subor- 
dinator whose Laplace exponent is given by 

p (A) = x A_(0(A + p) + g). (2.6) 

Moreover, the density h solves the equation 

POO 

(l-cx)h(x)= Tl p (\ogy/x)h(y)dy xG(0,l/c), (2.7) 



where IT p (z) = IL(z)e pz + qe pz . 

We remark that the transformation studied in Chazal et al. [9J is more general than 



the one presented in (2.6) and that they applied such transformation to Levy processes 
with one-sided jumps. We also remark that the subordinator whose Laplace exponent is 
given by (f) p has an infinite lifetime in any case. 

Our next goal is to study the behavior of the density of I e near 0. When q = 0, we 
work with the following assumption: 

(A) The Levy measure U belongs to the class C a for some a > 0, that is to say that the 
tail Levy measure IT satisfies 

lim ^ + ^ = e~ ay , for all y e R. (2.8) 

x^oo U(x) 

Observe that regularly varying and subexponential tail Levy measures satisfy this as- 
sumption with a = and that convolution equivalent Levy measures are examples of 



Levy measures satisfying (2.8) for some index a > 



Theorem 2.6. Let q > and £ = —(, where ( is a subordinator such that when q = 
the Levy measure II satisfies assumption (A). The following asymptotic behaviour holds 
for the density function k of the exponential functional I eq . 

i) If q > 0, then 

k(x) — >■ q as x 10. 

ii) If q = 0, then E[/~ Q ] < oo and 

k(x) ~E[r a ]Tl(logl/x) as xiO. 

In the sequel we will assume that q = 0. The above result will help us to describe the 
behaviour at oo of the density of the exponential functional of a particular spectrally neg- 
ative Levy processes associated to the subordinator (. In order to explain such relation, we 

6 



need the following assumptions. Assume that U, the renewal measure of the subordinator 
(, is absolutely continuous with respect to the Lebesgue measure with density u which is 
non-increasing and convex. We also suppose that E[£i] < oo. According to Theorem 2 in 
Kyprianou and Rivero [13] there exists a spectrally negative Levy process Y = (Y t ,t > 0) 
that drifts to +oo, whose Laplace exponent is described by 

V(A) = A0*(A) = -£-,, for A > 0, 

0(A) 

where <fi* is the Laplace exponent of another subordinator and satisfies 

J(0,oo) V J 

where 

a*=(c+[ xlWf 1 c* = {° c>0,orn(0,oo) = oo, 

V 7(0,00) 7 ' I l/n(0,oo), c = 0andn(0,oo)<oo, 

and the Levy measure IT satisfies 

u(x) = c*l{ x= oy + q* + LT (x), for x > 0. 

Let 7^, be the exponential functional associated to —Y, i.e. 



oo 



/,,, = / e Yfl ds 



t/, = I e us, 

o 



and denote its density by fe^,. From the proof of Proposition 4 in Rivero [21] the density 



k$ satisfies 



kJx) = q*-k( - ) , for x>0. (2.9) 

x \x 



The following corollary give us the asymptotic behaviour at oo of the density of the 
exponential functional of —Y. 

Corollary 2.7. Suppose that £ is a subordinator satisfying assumption (A) and such 
that its renewal measure has a density which is non-increasing and convex and let Y be its 
associated spectrally negative Levy process defined as above. Then the following asymptotic 
behaviour holds for the density function kf, 

h (x) ~ q*E \r a ] -n(log x) as x -> oo . 
x 



3 Proofs 



Proof of Theorem \2.1\ We start the proof by showing that the function 

1 



h(t,x) :- 



X t 



-{t<T } 



is such that 



h(t,x) d£ = 1, 



t > 0, x> 



for x > 0. 



(3.1) 



Then the result follows from the identity (3.1) and the fact that 

h(t + s) = qQ l [h(s,X t )l {t<To} j, for s,t>0, 
which is a straightforward consequence of the Markov property. 



Let us prove (3.1). From the definition of X and the change of variables u = r(t/x) 
which implies that du = x _1 exp{— £L t /x \}dt, we get 

/ h(t,x)dt = q dtE x _1 exp | - C(t/x)\ 1 Mt/x)<oo} 

=qE I x ' iexp {-^ t/x) y^ xfie ,i ds ft 

' POO 

= qE / l {um du =qE(P) = l. 



Uo 



We now prove that 



/ h(s)ds = F(l eq >t\ t>0. 



Indeed, let t > making a change of variables, using the semi-group property, and Fubini's 
theorem we have 



/oo y*oo / /»oo \ 

h(s)ds= h(s + t,l)ds = Q 1 ( h{s,X t )ds)l 

The result follows from the identity Qi(t < To) = P( I eq > t 



-{t<T } 



h(t<T ). 



a 



Proof of Corollary 2J2, Here, we use the same notation as above and we follow similar 
arguments as in the proofs of Lemma 5 and Proposition 1 in [3] . We first prove that for 
every < t < T and p > 0, the variable 

/■To l 
X? / -4vrds 



P+l 



a; 



is independent from a{X s , < s < t} and is distributed as 



~ p ^ds. 



As a consequence of the Markov property at time t, we only need to show that under Q x , 
the variable 



i" 



To 



o X. 



■p+i 



ds 



t ft 



is distributed as J^ q e p ^ 3 ds. Then the change of variables t = r(s/x), s = x J Q e^ u du 
yields 



xs 



T„ 



X, 



P+l 



T 



1 f t 

—ds = x' 1 / e ~ (p+1) ^w*)ds 



e -(p+D« t T e C t T dt 



- p «'dt, 



which implies the desired identity in law since (£j,0 < t < 0) and (^t,0 < t < e q ) have 
the same law. Hence, we have 



To 



xr 1 



ds 



X- p ;t<T 



<j>{p) + q 



which implies 



xr p ;t<T 



dt 



"(0(P) + 



Xr (p+1) ;t<T 



By iteration, we get that the function t H- Qi X t p ;t<To is completely monotone and 
takes value 1 for t = 0. Thus taking p = 1, we deduce that h(t) is completely monotone 
on (0, oo) and that lim^o h(t) = q. Finally from Theorem 51.6 and Proposition 51.8 in 
|23j . we have that the law of I e is a mixture of exponentials. □ 

Proof of Proposition 2. 3. The proof follows from the identity 

n 



E 



I? 



<p{n) + q 
Indeed, on the one hand it is clear that 



-E 



rn- 1 



n > 



E 



OO POO 

x n k(x)dx — n dyy 

n Jo 



n-l 



k(x)dx. 



On the other hand from the identity (2.3) with f(x) = e nx and a change of variables, we 



get 



n 



in) + q 



-E 



jn— 1 



oo 



i-OO 

-nx i „,n—l 



n\ U q (dx)e~ nx y n - l k(y)dy 
o Jo 

POO 

n I U q {dx) / ^e-^fc^dj/ 



o 







oo 







oo 



n / UJdx) / ^A^Ods 



o 



'o Jo 



Then putting the pieces together, we have 



dyy 



n-l 



k(x)dx~- I dyy 11 - 1 k(ye~ x )U q (dx) 
'o Jo 



for n > 0, 



which implies the desired result because the density 

1 



y ^ 



E(/ € 



q J Jy 



k(x)dx, 



is determined by its entire moments, which in turn is an easy consequence of the fact that 
k is so. □ 



Proof of Theorem 2.Jf., By Theorem 2.1 (when q > 0) and Theorem 3.9 in [2] (when q = 0), 
we know that there exists a density of 7 e > for g > 0, that we denote by /i. Moreover, in 
[7] it has been proved that the moments of I eq are given by 



E 



??.! 



nr=i (?+^; 



neN 



(3.2) 



where the product is understood as 1 when n = 0. 

We first prove that the function h : (0, oo) — > (0, oo) defined via 



h(x) 



cxh(x) 
0. 



Tl(\og(y/x))h(y)dy + q \ h(y)dy, if x G (0, 1/c), 

elsewhere, 



is a density for the law of J e and hence that h = h a.e. Then we prove that the equality 



(|2.5|) holds. In order to do so, it is enough to verify that 

x n h(x)dx 



nl 



n G N, 



niLi [q + <k 

since the law of I e is determined by its entire moments. 



10 



Indeed, elementary computations, identity (2.4) and the fact that 

6>0, 



e~ 9y U q (dy) 



1 



+ q 



give that for any integer n > 0, 



x n h(x)dx = c dxx n+ h(x) + / dxx n / dyH(log(y/x))h(y) 



+ q / dxx 11 / h(xe y )UJdy) 



+ dyh{y) / dxx n U{\og{y/x)) 



K=l (? + </•«) 



n!(n + l)c 

iEi 1 (? + #o) 



7i!(n + l)c 



/*oo /*oo 

+ g / t/g(dy) / dxx"/i(xe 4 
Jo Jo 



dyh(y)y n+1 / d^ e - (n+1)2 n(z) 

n JO 



+ g / U q (dy)e- (n+1)y / d^ n /i(z) 



n 



(1 - e-^)ll(dz) 



eg 1 (? + #o) nri 1 (? + 



rc+l 



n! 



z ^o 



r? 



7?! 



EKU0 + 

/■oo 

+ l)c+ / (l-e- (n+1)z )n(d^) + g 



[/ g (d2/)e-(" +1 ^ 



g + 0(n + 1] 



nili [q + M 



Now, let jV = jxGR: h(x) ^ h(x) > . By the above arguments, we know that the 
Lebesgue measure of N is zero. Let k : (0, oo) — ¥ (0, oo) be the function defined by 



'h(x) 
k(x) = { 1 



1 — ex 



if x e Af c . 
U(log{y/x))h(y)dy + q I h(y)dy ) , if x eAf. 



We now prove that k(x) satisfies equation (2.5) everywhere. If x G M c then we have that 
k(x) = h(x) = h(x), and hence equation (2.5) is verified. Finally, if x E Af, we have the 



11 



following equalities 

t*co /»oo 

cxk(x) + / n(log(y / x))k(y)dy + q / k(y)dy 

J X J X 

/*oo poo 

= cxk(x) + U(\og(y/x))k(y)l {yeA fc } dy + q k(y)l {yeAf c } dy 

J X J X 

/*00 /*00 

= cxk(x)+j U(\og(y/x))h(y)l {yeA fc } dy + q h(y)l {yeAf c } dy 
U(hg(y/x))h(y)dy + q h(y)dy) ) 



ex 



1 — ex 



k(x). 



x 

oo 






+ 



U(\og(y/x))h{y)dy + q / h{y)dy 



Conversely if k is a density on (0, 1/c) satisfying equation (2.5), from the above compu- 



tations it is clear that k and I e have the same entire moments. This implies that the k 
is a density of the exponential functional I eq . □ 

Proof of Theorem 2.6. The proof consists of three steps. First we show that when q = 0, 
E[/~ a ] < oo, then for q > we obtain a technical estimate on the maximal growth of 
k{x) as x \. 0, and finally the statement of the theorem. 

Step 1. Here we assume that q = and prove that E[/~ a ] < oo. The case a = 
is obvious. For a G (0,1), we have from Theorem 2 in [1] that there exists a random 
variable R, independent of £, such that IR = e, where e follows a unit mean exponential 
distribution. Since E[e _a ] < oo, the result follows. 



Finally let a > 1. With (2.5) and some standard computations, we find 

> ^oo ^oo /*oo 

x^ l3 ^ 1 k(x) dx = c / dxx~ l3 k(x) + / dxx^ 13 ^ 1 I dyTl(\og(y/x))k(y) 



o 



cE 
cE 



J-/3 



+ 







dyk(y) / dsa; ^ 1 Il(log(?//x)) 



o 



/ dyy- p k{y) / due^IIfa) 
o Jo 



-Ie 



J-/3 



that is to say 



E 



/-/J-i 



-c/3+ / fi-e^jn(d^; 

0(-/3) 



E 



I'? 



■P 



(3.3) 



where <j> is the Laplace exponent of £, which can be extended to (—a, oo) since for f3 < a 



/OO /»O0 

(e^ - 1) n(du) =p n(log(^))^ 1 dz < oo. 



(3.4) 



12 



To see that (3.4) holds, note that n(log(z)) is regularly varying with index —a by (2.8). 
Hence II(log2) = z~ a £(z) for a slowly varying function C. and we can apply Proposition 
1.5.10 from Bingham et al. [3]. 



Now, by iteratively using (3.3) we see that for E[/~ a ] < oo it is enough to have 
E[/~ a ] < oo for some a 1 G [0, 1). But this obviously holds if a' = 0, while if oi G (0, 1) it 
holds by the same argument as we used above for the case a G (0, 1). 

Step 2. We assume that q > 0. For q = 0, let p be any function such that p(0) = and 
min{a — 1, 0} < p(a) < a for all a > 0. When q > the function p will be taken as zero 
and hence the symbol p(a) will be taken as 0. The goal of this step is to show 



k(x) 



x 



p(a) 



stays bounded as x \. 0. 



(3.5) 



k(x)/x p ( a \ We can write (2.5) as 



Observe that when q > it follows from (2.5) that liminf^o k(x) > q. Set h(x) : 

qx p( - a) F(I e > x) 



1 — ex 



X 



IHa ) h(xz) 

h(x) h(x) 



n(iog(z))^ 



dz + 



(3-6) 



We argue by contradiction. Take some x G (0, 1/c). If h were not bounded at 0+, then 
l{x<x}h(x) would keep on attaining new maxima as x 4- 0. (Note that x is present just 
to make sure this statement also holds if k is not bounded at 1/c—.) In particular this 
means that a sequence of points (a; n ) n >o exists with x n \. as n — > oo and such that 
h(x n ) > sxvp xe ! Xn ^ h(x). We will show that this implies 



- ' i. 



n(iog(z))^ Q )^^ 



h(x r 



dz + 



qx p n ia) F(I eq > x r 
h(x n ) 



^0 



as n 



oo, 



which indeed contradicts (3.6) since 1 — cx n — > 1 as n — > oo. Observe that if q > and 



h is not bounded at 0+ then the second term in the latter equation tends to 0, because 
p{ot) = 0, by construction. So we just have to prove that the first term in the latter 
equation tends to 0. For this, we have 



- ' r. 



n(iog(z))^ (a) ^^ 



h(x r . 



dz 






n(hg(z))z 



,,< a) h{x n Z) 

h(x n ) 



dz 



+ %r, 



n(iog(^))z p ^^^ 



(3.7) 



X/Xfi 



h{x n ) 



dz. 



We first deal with the first integral on the right hand side of (3.7). By construction of the 
sequence (x n ) n > , we have h{x n z) < h(x n ) for any z G [l,x/x n ], hence 



,./,.,■ 



•E / -En 



n(to g (*))**°>^r^ 

h{x n ) 



dz < x r 



■E/ -^n 



U(log(z))z p(a) dz. 



(3.J 
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If q > or a = (recall p(0) = 0), we can take any 1 < zq and write 



"£/£n rzo r&/x n 

.;•„ / U(\og(z))dz = x n / U(\og(z)) dz + x n / II(log(z)) dz 



< av 



in 



n(log(z))dz + x r 



X 



X i. 



z U(\og(z )), 



where the inequality uses that II is decreasing. Letting n — > oo, recalling that x n \, 0, we 
see that the first integral on the right hand side vanishes while the second term tends to 
xTMAogzo). As we can make this term arbitrarily small by choosing zq large enough, since 
n(logz) — > as z — > oo, it follows indeed that (3.8) vanishes. 

Next, let a > 0. Since a — 1 < p(a) < a, we can choose some f3 G (0, a) such that 
p(a) — (3 + 1 G (0, 1). Using this we find 



• < i 



/x/x„ rx/x n 

Tl(hg(z))z p{a) dz = x n Ii{\og(z))z p - 1 z p{a) - p+1 dz 



< X r 



X 



,1 ;■ 



p(a)-/3+l r x/x n 



Uilogiz^z^dz, 



and the right hand side indeed vanishes as n — > oo, again since x n I and by (3.4). 

It remains to show that the second integral on the right hand side of (3.7) vanishes as 
n — > oo. We have 

x n [^ U(\og( z ))z p ^ h p^dz<x n U(log(x/x n ))-^— r z^h(x n z)dz 

Jx/xn "'(En) il{%n) Jx/x n 

U{\og{x/x n )) 1 ' v 



X P n {a) h( Xn ) 



k{u) du, 



where the inequality uses that II is decreasing and the equality the definition of h together 
with the substitution u = x n z. Since k is a density and by assumption h(x n ) — > oo as n 
goes to oo, for the right hand side to vanish it remains to show that Ii{\ogx/x n )/xn 
stays bounded as n increases. If q > or a = (recall p(0) = 0) it is immediate since IT 
is decreasing. If a > 0, for any 1 < zq < z integration by parts yields 



U{log{z))z pia) = p{a) / n(log(M))M p(a) " 1 du + 
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u 



p(a) 



p(a) 



dU(\og(u)) + U(\og(z ))z PW . 



zit 



Now, if we let z goes to oo, then since p(a) < a we see from (3.4) that the first integral 
in the right hand side stays bounded while the second integral is negative on account of 
the fact that IT is decreasing. Consequently the left hand side has to stay bounded and 
we are done. 



Step 3, case q = 0. Denote C a = E[J a ], which is finite by Step 1. From (2.5) we 
obtain for all x > 0, 



'1 — ex) 



k(x) 



Yi(\og(y/x)) 



n(io g (i/x)) J x n(io g (i/x)) 
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k(y) dy. 



(3.9) 



Using this equation together with k > 0, Fatou's lemma and the fact that II has an 

n(log(y/x)) 



exponential tail (cf. (2.8)) yields 
k(x) 



lim inf : 



lim inf 



cxk(x) 



+ lim inf 



*4o n(iog(i/x)) <4° n(iog(i/x)) ^o j x n(iog(i/x)) 



-k(y)dy 



> I y~ a k(y)dy = C a . 
'o 



On the other hand, for any e>0we have as x \. 



Jf 



n(log(y/a;)) 
n(log(l/x))' 



%)d|/-> / y- a %)dy<C a . 



If a > 0, this follows from the fact that the convergence (2.8) is uniform over y e [e, oo), 
see e.g. Theorem 1.5.2 in [5]. If a = this uniformity holds only over intervals of the form 
[e, Xq], in which case we can write the left hand side as the sum of integrals over [e, Xo] 
and [xq, oo), the former in the limit again is bounded above by C a , while for the latter 
we can use that IT is decreasing to see 



U{log{y/x)) n(log(zo/aO) 

-k{y) dy < 



xo n(io g (i/s)) 



n(io g (i/x)) 



Kv) d v> 



then letting first x — >■ oo, thereby using (2.8), and then xq — > oo it follows that this term 
vanishes. 

So it remains to show that 



lim sup 



U(log(y/x)) 



n(io g (i/x)) 

For this, we get for e small enough and x < e: 



k(y) dy — > as e — > 0. 



n(io g (i/x)) 



U(log(y/x))k(y)dy 



x 



< 



n(iog(i/x)) 

Cx 

n(b g (i/x)) 

n(io g (i/x)) 

CV+p(a) 



e/i 



n(log(^))fc(a;^) d^ 



e/rc 



e/a 



n(iog(z))(^) p(Q) d^ 
n(iog(^))z p(a) d2 



n(log(l/x)) Vx7 



£\ p(o) + l 

i Il(log(e/x)) as x J, 0, 



where C and C" are constants, the inequality holds by Step 2 (cf. (3.5)) and the asymp- 



totics follows from Karamata's theorem (see e.g. Theorem 1.5.11 in [5]), which indeed 



applies here since n(log(z)) is regularly varying with index —a (cf. (2.8)) and by con- 



struction (see Step 2) p(a) > a — 1. Now, using (2.8) we see that the ultimate right hand 
side goes to C'e p( - a ^ +1 ~ a as x J, and this indeed vanishes as e — > since by construction 
p(a) + 1 - a > for all a > 0. 



15 



Step 3, case q > 0. We will prove that 

U(log(y/x))k(y)dy — -> 0. 

x— >0 

By Step 2, we can assume that k is bounded by K > q, in a neighborhood of + . Let 
5 > 1 fixed, for x small enough we have that 

x8 pxS 

U{log{y/x))k{y)dy <K U(log(y/x))dy 



/•log S _ 

= K Il(u)xe u du 

Jo 

/•log <5 _ 

< KxS / U(u)du y 0. 

Jo ^° 



Also, we have that 



U{log{y/x))k{y)dy < U{hgS) / k(y)dy — f U(\ogS) 



xS JxS x ^° 



We conclude by making § — > oo. Indeed, using equation (2.5) and the above arguments 
we conclude that 

(1 - cx)k(x) - q¥(I eq > x) ► 0, 

x— >0 

from where the result follows. 

□ 



4 Examples and some numerics 



In this section, we illustrate Theorem 2.3 , Corollary 2.4 and equation (2.9) with some 
examples and we also provide some applications of Theorem 2.5. 

Example 1. Let q > and consider the case when the subordinator is just a linear 
drift with c > 0. By a simple Laplace inversion, we deduce u q {x) = c~ 1 e~^ x . Thus, from 



identities (2.5) and (2.4) we get 



(1 - cx)k{x) = - k(xe y )e~i y dy, x e (0, 1/c). 

C J[0,oo) 

After straightforward computations, we deduce that the density of I e is of the form 

k(x) = q(l-cx)c~ 1 , xe(0,l/c). 

Let p > and note that 

6 , q p(p + l)T(q/c) 



,(9) = c6 + q- and c 



py ' H 9 + p ' p c^Tip + q/c+l) 
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According to Corollary 2.4, the density of the exponential functional of the subordinator 
whose Laplace exponent is given by P , satisfies 

h(x) = ^ r^+nrfaAO ^ 1 "" cx)q ' c ~ 1 ' for x e (0, 1/c) ' 

in other words the exponential functional has the same law as c~ 1 B(p + l,q/c), where 
B(p + 1, q/c) is a beta r. v. with parameters (p + 1, q/c). 

Now, let us consider the associated spectrally negative Levy process Y whose Laplace 
exponent is written as follows 

«,(A)- y - A < A+ "> 



,(A) c(\ + p) + q 



From (2.9), we deduce that the density of the exponential functional /</, associated to Y 
satisfies 

cp + qT(p + l)T(q/c) 
Hence 1^ has the same law as c(B(p, q/c))^ 1 . 

Example 2. Let q = c = 0, (3 > and 

where a G (0, 1] and s > a. Thus, the Laplace exponent has the form 

0T(a(e - 1) + s) 



<K0) = p- 



T(a9 + s) 



In this case, the equation (2.5) can be written as follows 



f3 f°° _2-i / i \a~l 

fe(^) = r(a + 1 x / {y/x) a [(y/x) a - ij Hy)dy 



fix 
Via) 



[ (z + iy-'z^kfxiz + l^dz, 



where we are using the change of variable z = (y/x)* — 1. After some computations we 
deduce that 

k( z ) = J^- z ^ e -^ z) \ for z>0. (4.10) 

aT(s) 

In other words I has the same law as /3 _1 7", where 7 S is a gamma r.v. with parameter s. 
If a = 1, the process £ is a compound Poisson process of parameter (3 > with 



exponential jumps of mean (s — 1) x > 0. From (4.10), it is clear that the law of its 
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associated exponential functional has the same law as 7( s „8), a gamma r.v. with parameters 

{8,0). 

We now consider the associated spectrally negative Levy process Y whose Laplace 

exponent satisfies 

_A*_ Ar(aA + s) 

^ >~ 0(A) ~ /3T(a(\-l) + s)' 
The density of the exponential functional 1$ associated to Y is given by 

h{x)= -f^ x -" a e-W°> 1/a , x>0. 
ai (s — a) 

We remark that when a = 1, the process £ is a Brownian motion with drift and that 
the exponential functional 1^ has the same law as r y7 s 1 _ 1 ay This identity in law has been 
established by Dufresne [TO] . 
Next, let p > and note that 

U0)=^ 9 ( : P \ l) \ S) and ~-^P + s) 



T(a(6 + p) + s) p pPT(s) ' 

According to Corollary 2.4, the density of the exponential functional of the subordinator 
whose Laplace exponent is given by (f) p , satisfies 

a{s+ap)/a 
h{ X ) = —, (ap+s~a)/a e -(M I- for % > Q 

aT(ap + s) 

i.e. it has the same law as /3 _1 7" +s . In particular, the density of the exponential functional 
of its associated spectrally negative Levy process satisfies 

o(s+ap-a)/a 
k,(r)- — — T -(ap+s)/a -{(3/xY^ n 

aT{a{p — 1) + s) 
Example 3. Finally, let a 6 (0, 1), /9 > a, c = 0, q = T(0)/T(p - a), 

— 1 f 00 P (l+a-P)x/a 1 , x a -i 

II(z) = — - — - / -^- -—da; and uJz) = s e~^~ 1)z/a ( e« - 1 ) . 

v; T(l-a)J z (e x / a -iy+ a qK ' r(o+l) V ) 

The process £ with such characteristics is a killed Lamperti subordinator with parameters 
(l/r(l — a), l + a — /3, 1/a, a), see Section 3.2 in Kuznetsov et al. [15] for a proper definition. 
From Theorem 1.3 the density of I e satisfies the equation 

k(x) = / w r / T^ tt— dx + ^F^ r ( e- - 1 J fcfxe^dy. 

w 7 \T(l-a)J z ( e */a_i)i+« r(/3-a)r(a + l)V 7 ^ l ;y 

Since the above equation seems difficult to solve, we use the method of moments in order 
to determine the law of I e . We first note that 



E 



4 n 

e <3 



n\r(p) 

' T(an + (3) 
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and that in the case /3 = 1, the exponential functional I eq has the same distribution as 
X~ a , where X a is a stable random variable, i.e. 



E 



-\x a 



exp{-A a }, A > 0, 



see Section 3 in [3J. Recall that the negative moments of X a are given by 

r(l + n/a) 



E 



x: n 



r(i + n) 



n > 0. 



Now we introduce Lum and A, two independent r.v. whose laws are described as follows, 



P(L M) G dy) = E 



oro?) i 
r(p/ a )xtxi 



edy 



and 



P(A G dy) = (p/a - l)(l - xf /a - 2 l m {x)dx. 



It is important to note from example 1, that A has the same law as the exponential 
functional associated to the subordinator a which is defined as follows 

<r t = t + P/a-\, t > 0. 

On the one hand, it is clear that 



E 



>,/3) 



aT{p) 



E 



v-(an+/3) 



r(/3) Tjn + P/a) 
T(P/a) T(an + (3) 



and on the other hand, we have 



E 



.4' 



r(n + l)T(p/a) 



T(n + /3/a) ' 
which implies the I eq has the same law as L a ^A. 

Finally we numerically illustrate the density k and its asymptotic behaviour at for 
some particular subordinators (. Let us first shortly discuss the method we used. Clearly 



the equation (1.2) motivates the following straightforward discretisation procedure: ap- 



proximate A; by a step function k, i.e. 



k(x) 



N-l 

/ j ^-{xe[xi,x i+ i)}Vi} 



i=0 



where = Xq < x\ < . . . < x^ = 1/c forms a grid on the x-axis. The heights yi can then 



be found by iterating over i — N — 1, . . . , 0, thereby at each step using (1.2) with x = Xi 
and k replaced by k. Two remarks are in place here. 
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Firstly, as (1.2) is linear in k the condition that fc is a density is required to uniquely 
determine the solution. This translates to the fact that the numerical procedure discussed 
above requires a starting point, i.e. the value yjv-i > should be known. (Of course, 
starting with y N = yields k = 0.) We proceed by first leaving y^-i undetermined, run 
the iteration so that every j/j in fact becomes a linear function of i/n-i, and then find Dn-i 
by requiring that k integrates to 1. 

The second remark is that even though any choice of grid would in principle work, we 
found one in particular to be useful. Indeed, if we set x n = (l/c)A N ~ n for some A less 



than (but typically very close to) 1, equation (2.5) yields the following relation 



cc 



1 - cx n )y n = / U(\og(y/x n ))k(y) dy = x n U(\og(z))k(x n z) dz 

J x n J \ 

N-l „oo N-l A n-i-l 

n Y/£/i / n(\og(z))l {XnZe{x . tXi+l)} vz = x n s ^y i / n(log(z)) dz. 

TZ Ji TZ M"-' 



i=n 



The approximation this setup yields is very efficient in comparison with e.g. the approx- 
imation using a standard equidistant grid, due to the fact that in this case we need to 
evaluate only iV different integrals numerically 

First we consider two examples for which the density k of / is explicitly known. The 



first one is taken from Example 2 with a — 1, f3 — 2 and s = 3/2. In this case from (4.10) 
we have 

2 5/2 

k(x) = —=x 1/2 e 2x for x > 0. 
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See Figures 1-4 for plots of the density k, the difference k — k (where k is obtained by 
the above method with A = 0.998, yielding a grid of ps 4500 points and a few min- 
utes computation time on an average laptop), the ratio fc(x)/II(log(l/x)) and the ratio 
k(x)/U(\og(l/x)) respectively. 

The second explicit example is taken again from Example 2 with /3 = 1 and s = 1 and 
a = 1/2. In this case from (4.10[), we have 
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k(x) = 2xe x for x > 0. 

It is important to note that II satisfies (A) with a = 1. In this case, Figures 5 -8 show 
plots plots of the density k, the difference k — k, the ratio k(x) /Il(log(l / 'x)) and the ratio 
A;(a;)/n(log(l/a;)) respectively. 

Next we look at two examples where no formula for k is available. The first one is 
when £ is a stable subordinator with drift, i.e. c = 1 and n(dx) = x~ 1 ~ a dx, where we 
take a = 1/4. See Figures 9 & 10 for a plot of k and the ratio A;(x)/II(log(l/x)) respectively. 



Note that this is an example of a Levy measure satisfying (2.8) with parameter 0. 

Finally, the second example is a subordinator £ with zero drift and Levy measure of 
the form Il(dx) = a; _1//4 exp(— x n ) dx. Figure 11 shows k for n = 1 (blue), n = 2 (purple) 

* All computations were done in the open source computer algebra system SAGE: www.sagematli.org 
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and n = 3 (green) respectively. Figure 12 shows the ratio k(x)/U(logl/x) for the case 



n—1, since then (A) is satisfied with a — 1. 
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Figure 2: The difference k — k 




Figure 3: The ratio k(x)/U()og 1/x) 



Figure 4: The ratio A;(a;)/n(logl/x) 





Figure 5: The density function k 



Figure 6: The difference k — k 
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Figure 7: The ratio k(x)/U()ogl/x) 



Figure 8: The ratio k(x)/H(logl/x) 





Figure 10: The ratio fe(x)/II(log 1/x) 



Figure 11: Density functions k 



Figure 12: The ratio k(x)/U(logl/x) 
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